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ABSTRACT. 

We present the results of the first two color VR observations of the Algol- 
type binary system WZ Crv (12 fc 44 m 15.19", -21 d 25 m 35.4 s (2000)) which were 
obtained using the remotely controlled telescope TOA-150 of Tzec Maun Ob- 
servatory. We determined the moments of individual minima, the orbital pe- 
riod and its derivative, the initial epoch, color indices V-R and temperature 
estimates of the components. Also we noticed that the phase curve is asym- 
metric: the second maximum is higher than the first one. It indicates that 
there is a spot in the photosphere of one of the stars in this system. 

Keywords: Stars, Eclipsing binaries, Algol-type, individual: WZ Crv 
PACS: 97.80.Hn Eclipsing binaries; 97.10.Gz Accretion and accretion disks 
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I INTRODUCTION. 



The Algol-type binary system WZ Crv had been discovered in 1936 and published 
by Luyten [lj in 1937 with the number 686. Then it got the preliminary numbers 
NSV 05914, EC 12416-2109 and GSC 6110.00930. The primary minimum is rather 
deep. According to Luyten the brightness range is from 13. m to 15. m 0. In the 
General Catalog of Variable Stars (Samus' et al. [2J) the value of the period of WZ 
Crv is P = l. d 78878, and the range of the brightness variations is 12. m 76 - 14. m 7 in 
V-band. During the Edinburg-Cape Blue Object Survey (Kilkenny et al. |3j) the 
spectral classes G and F7 for the components had been determined. Two minima 
timings were observed by Otero & Dubovsky [4J (CCD observations) and Locher [5] 
(visual observations). 



II OBSERVATIONS 

We observed the south Algol-type binary system WZ Crv (12 h 44 m 15.19 s , 
-21 d 25 m 35.4 s ) = NSV 5914 during the period from JD 2455321 to JD 2455395 using 
the remotely controlled astrophysical refractor TOA-150 of Tzec Maun Observatory 
(D = 150mm, F = 1100mm), located in Moorook, South Australia. This telescope 
was equipped with SBIG ST-L-6K CCD camera with a filter wheel. The field of 
view is 57.5' x 86.3', the pixel scale is 3.37" /pixel. We observed in photometric 
V-Bessel and R-Bessel filters and got 317 fit points in each filter. The 120 sec 
exposures were used. The images were calibrated by dark frames and flat field 
frames. 

To test our photometric system we observed the field of the variable star AE 
Aqr. The BVRI magnitudes in standard bands of 228 stars in the vicinity of AE 
Aqr were determined by Henden [5]. We used 93 most bright stars and determined 
the color transformation coefficients: 

d(Ri - R) 

- 0.015 ±0.035; 



d(V - R) 
d(Vj - V) 
d(V - R) 



-0.033 ±0.043. 



where Ri and Vi are the magnitudes in the instrumental system, and V and R are 
the standard magnitudes (by Henden). This means that our photometry system is 
close to the standard one. 

There is no photometry data for reference stars in the field of WZ Crv in 
AAVSO database or in any previous papers about this star. But there are SDSS 
[7] photometry data in the "u, g, r, i, z" photometric system for the stars in the 
vicinity of WZ Crv. We used 3 reference stars to calibrate our photometry: USNO- 
A2.0 0675-12226641 (12 h 44 m 54.624 s , -21 d 30 m 27.21 s ), USNO-A2.0 0675-12232415 
(12 h 45 m 35.274 s , -21 d 33 m 42.93 s ), USNO-A2.0 0675-12233718 (12 h 45 m 44.211 s , - 
21 d 16 m 07.06 s ). The transformation formulae (Lupton, 2005, 0) were used to 
convert the "u, g, r, i, z" data into BVR system. The information about the 
reference stars is given in the Table [TJ Their positions and the position of the 
investigated star WZ Crv are shown on the Fig. [TJ 
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Table 1. 



# 


USNO-A2.0 


u 


g 


r 


i 


z 


B 


V 


R 


1 


0675-12226641 


15.219 


13.870 


13.295 


13.230 


13.340 


14.277 


13.534 


13.132 


2 


0675-12232415 


15.781 


13.953 


13.248 


13.504 


12.818 


14.401 


13.541 


13.179 


3 


0675-12233718 


16.534 


14.369 


13.532 


13.284 


13.046 


14.858 


13.881 


13.315 



N 
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Figure 1. Finding chart. The field of view is 30'x25'. 


The positions if three reference 


stars and WZ Crv are marked. North is up, East is 


left. 


III Data analyzing. 





We determined the period, the initial epoch and the values of maxima and minima 
in both, V and R, bands with corresponding errors using our observations. The 
software FDCN (Andronov [HI [10]) was used to compute the coefficients of the 
statistically optimal trigonometric polynomial fit, using the least squares method, 
which allows to calculate the necessary parameters and their errors estimates. The 
statistically optimal degree of the trigonometric polynomial s = 9. The determined 
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values are: 



P = l d 788788 ± d 000033 

T = HJD2455354.63500 ± 0.00031 



min / ( J R) = 14 m 004 ± m 005 
mm n (R) = 12 m 751 ± m 006 
Max/(i2) = 12 m 588 ± m 005 
Max J/ ( J R) = 12 m 534 ± m 005 



miri/(y) 
mm a (V) 
Max/(y) 
Max /7 (V) 



14 m 564 ± m 009 



13 m 052 ± 0. m 010 



12 m 889 ± m 006 



12 m 868 ± m 007. 



The phase curves are shown on the Fig. [2j Besides, we made the V-R phase 
curve; the result is presented on the same Figure in the same scale. It shows that in 
the primary minimum the temperature is lowest, and in the phase of the secondary 
minimum the temperature is highest. 

IV Analyzing of the Phase Curves 

At the main (deeper) minimum, there is a horizontal part which is typically inter- 
preted as a total eclipse of a smaller star (called "primary") by a larger one (called 
"secondary"), so during this time interval we see only a larger star. The brightness of 
the star is m 2 = m min i. To determine this value, we have computed a sample mean 
value of the measurements during this phase interval of assumed constant brightness 
either in the filter V, or R. To determine the brightness of the eclipsed star, we 
have to use the expression Eq = E\ + E2, where Eq, E\ and E2 are intensities of the 
binary system (both stars), primary and secondary stars, respectively. The intensi- 
ties E are related to brightness (stellar magnitudes) m using the classical Pogson's 
formulae: 



where E$ and ms are the intensity and brightness of some calibration source. From 
these expressions, one may determine the magnitude of the eclipsed primary star as 



The value of mo is a subject of a special discussion. Sometimes the out-of-eclipse 
parts of the light curve are not constant due to the "reflection" effect (heating of 
atmosphere of one star by an emission of another star) and the "ellipticity" effect 
(due to a tidal distortion of one or both stars). From the mathematical point of 
view, these effects may be modeled as a second-order trigonometric polynomial 

x(0) = d + C 2 cos(2tt0) + C 3 cos(4vr0) + C 4 sin(2vr0) + C 5 sin(4vr0). 

The coefficients C4 and C5 are close to zero for symmetric light curves. However, 
for WZ Crv, the maximum at phase 0.75 is brighter than that at the opposite phase 




m-ms = -2.5 lg — , 



m\ = mo — 2.5 lg (l — 10 



0.4(m2-mo) 
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0.25, so C4 is distinctly different from zero, whereas C5 = within error estimates. 
Thus we decreased the number of the determined parameters to 4. The formulae 
for the error estimates were presented by e.g. Andronov [9j HD]- The shape of 
the secondary minimum shows a linear-like descending and ascending branch linked 
with a parabola. Thus we have used the "asymptotic parabola" fit (Marsakova & 
Andronov Results for the brightness maxima in both filters V, R are presented 

above. We estimated the smoothed brightness m at the phase = of V = 
12 m 953 ± m 007 and R = 12 m 634 ± m 004. For <p = 0.5, V = 12 m 917 ± m 006 and 
R = 12 m 601 ± m 004. Thus the brightness of the primary is V x = 13 m 233 ± m 014, 
Ri = 13 m 000±0. m 008 and (V-R)i = m 237±0 m 016. Similarly, we determined the 
brightness of the eclipsed part of the secondary corresponding to the phase (f> = 0.5: 
V 2e = 15 m 251 ± m 014, R 2e = 14 m 824 ± m 054 and (V - R) 2e = m 427 ± m 109. 
The ratio of the intensities from the secondary to the eclipsing part is E 2 /E 2e = 
1.94 ±0.20 and E 2 /E 2e = 2.19 ±0.05 for V and R. These values are the same within 
their error estimates. If interpreting this ratio as the ratio of the eclipsed surfaces 
(in the simplest model of zero limb darkening), one may estimate the ratio of the 
radii R 2 /Ri = (E 2 /E 2e ) 1 / 2 of 1.40 ± 0.07 and 1.48 ± 0.04, respectively. Theoretical 
achievements are that this ratio is constant. However, the slight difference between 
the ratio E 2 /E 2e for V and R may be interpreted as a larger coefficient of the limb 
darkening in V as compared with that in R. This is in an agreement with results of 
Al-Naimiy [J2]. 

Using the color indices V — R, we estimated the temperatures of the primary and 
secondary stars in this binary system. In the primary minimum the cooler bigger star 
eclipses the hotter star. The flatness of the primary minimum (Fig. [3]) indicates that 
in the primary minimum the eclipse is total. Thus, it is possible to determine the 
approximate temperature of the first component, using the color index V — R in the 
phase of primary minimum. We determined the corresponding temperature, using 
the empirical formula (Andronov & Baklanov |13j). which describes the dependence 
between the color index V — R and the spectral class. There is a misprint in the 
English version of the paper, the correct version is: 

= 1.948(±0.011) + 1.474(±0.019) ■ (V - R - 0.681). 

In the total primary eclipse the value of the color index is V — R = 0. m 560 ± 
0. m 012. Thus, the temperature of the second component is T = (5650 ± 66)7^. 
Using the calibration of spectral class temperatures of Allen (1973) |14] the "3-ct" 
temperature "corridor" corresponds to the range of spectral classes of G2 - G6. Using 
the determined color index (V — R)i = 0. m 237 ± 0. m 016 for the primary component, 
we estimated the temperature of first component (7874±166)-fT and range of spectral 
classes A6 - Fl. 

Assuming that the inclination of the orbit is rather close to 90°, we determined 
the lower limits for the relative radiuses of the stars, using the formulae 

Ri sin(irD) — sin(7r<f) 
~a ~ 2 

R 2 sm(nD) ± sin(7rd) 
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where R\ and R 2 are the average radiuses of the stars, a is the orbital separation, i. 
e. a distance between the centers of the components, d is the duration of the total 
eclipse (in phases), D is the total duration of eclipse, from the beginning till the 
very end (in phases too). The phases of the eclipse are from —0.0915 to 0.0915, the 
phases of the total eclipse are from —0.017 to 0.017, which means that the values of 
D and d are: 

D = 0.183 
d = 0.034 

Thus, the lower limits of the relative radiuses: 

Ri 

— = 0.219 

a 

R2 

— = 0.325 

a 

The ratio R 2 /R\ = 1.48 determined from geometrical constrains is an excellent 
agreement with that obtained using the fluxes. Thus we conclude that the inclination 
angle is close to 90°. 

The most interesting feature of the phase curves is the statistically significant 
difference between values of maxima in both filters. This difference reaches 0. m 060 ± 
m 003 in filter R and m 039 ± m 005 in filter V. The cause of such effect might be 
a cold or a hot spot in the photosphere of one of the stars. 

To find out if the spot is steady, we checked the ASAS database. There are 
493 fit points in filter V, which were obtained during the period from JD2451876 to 
JD2455049. Using these observations we determined the period and the initial epoch. 
The values are P = l. d 788797 ± d 0000021 and T = #JL>2453494.2755 ± 0.0014 
respectively. The phase curve is shown on the Fig. 0J There is no statistically sig- 
nificant asymmetry on this curve, perhaps, because of rather big error of individual 
observations. Thus, using ASAS data, we can't conclude on the stability of the spot. 

On the other hand, according to the previous calculations, the spectral classes 
of the components of WZ Crv are G2 - G6 and A6 - Fl, which nearly completely 
excludes the possibility of forming rather steady big cold spots in the photosphere 
of at least one of the stars. More possible solution is the presenting of a hot spot as 
the result of deep impact of the accretion stream from the first component (donor) 
directly into the photosphere of the second one (accretor) without forming the accre- 
tion disk. In this case the hot spot will be rather steady, which might be confirmed 
by further observations. 

V Period Variations and Mass Transfer 

To study stability of the orbital period using the O-C analysis [15J, we compiled a list 
of the minima timings. The individual times of minima were published by Otero & 
Dubovsky (2004) and Locher (2005). added the median time of minima, calculated 
from our observations, from the NSVS database and from the ASAS database. ASAS 
data covers rather long period of observations (9 years), so we divided this array 
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into three parts (JD 2451876-2452879, JD 2452977-2453908, JD 2454091-2455049) 
and analyzed them individually. For the ASAS data, we have used several methods. 
The period and the initial epoch of the best "complete" trigonometric polynomial 
fit of order 9 were used for further analysis. As for these data there is no evidence 
for a statistical difference of maxima, we have used a "symmetric" trigonometric 
polynomial fit 

s 

x{4>) = d + Y^ C j+1 cos(2ttj(0 -<M), 

3=1 

where is the phase corresponding to the moment of time t: <j> = (t — Tq)/P — 
INT((t — Tq)/P), where To is an "initial epoch", P is a period and INT(x) is an 
integer part of x. 

The coefficients Cj were determined using a linear method of least squares 
and the parameter O (phase of the main minimum) was determined by itera- 
tions (non-linear method of least squares). Using the value of (p , we have com- 
puted a "mean" moments of time of minimum, which are closest to the mean time 
of the observations similar to a "complete" trigonometric polynomial fit x(<f>) = 
C\ + ^2j =1 Cj + i cos(27rj(0 — (f>oj)) with phases Oj -, which are individual for each har- 
monic. Symmetric trigonometric polynomial of 9-th order were used for the NSVS 
data and the 3 intervals of the ASAS. They are listed in the Table [2J 
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Figure 2. The phase curves of WZ Crv in R and V filters and the phase curve of 
the color index. 
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Figure 3. The shape of the primary minimum in the phase of full eclipse in R and 
V filters, obtained on JD 2455360 (2010.06.12). 
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WZCrv-ASAS data 
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Figure 4. WZ Crv phase curve according to ASAS data in V-band. 



Table 2. Minima timings. 



HJD 


HJD error 


E 


O-C 


source 


2451596.37218 


0.00103 








NSVS 


2452086.49800 


0.00800 


274 


-0.00209 


Otero 


2453440.59300 




1031 


-0.01961 


Locher 


2452378.07042 


0.00176 


437 


-0.00212 


ASAS 1 


2453483.53860 


0.00201 


1055 


-0.00492 


ASAS 2 


2454535.36234 


0.00242 


1643 


0.01148 


ASAS 3 


2455354.63500 


0.00031 


2101 


0.01923 


Virnina 



For all ASAS data, we obtained HJD min = 2453494.27648 ± 0.00115. However, 
this value corresponds to a wide interval and thus shifted from the smoothing curve. 
The minimum of Locher deviates significantly from the fit as well. This is due to 
a large error estimate, which is by a factor of 3-27 times larger than the minima 
timings determined from our observations and that of the NSVS and ASAS. Thus 
we excluded the values from the ASAS (three subintervals together) and of Locher 
[5] from our analysis. As there is no error estimate of the timing was published by 
Otero & Dubovsky [I], we have used the same weight for all minima timings. 

Thus, we analyzed 7 times of minima, and made the O-C diagram as the depen- 
dence between the cycle number and the residual O — C = Tg — (T + PE), i.e. the 
difference between the observed and calculated times of minima. The O-C diagram 
exhibits a parabolic shape, which is characteristic for secular period variations (cf. 
Tsessevich [16J). The smoothing parabola is shown on the Fig. [5j The least-squared 
method yields the equation for time of minima: 

T min = 2451596.3722+ 1.788795 -(E - 918) + 1.09 • 10~ 8 -(E - 918) 2 

±0.002 ±0.000002 ±0.3M0" 8 
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From the diagram and the equation it is clear that the orbital period and, accord- 
ing to Kepler's law, the distance between the components increases. In assumption 
of conservative system, we determined dP/dt: 

dP/dt = (1.2 ± 0.3) • 10~ 8 days/yr, 

and the characteristic time of the period variations is 

P/(dP/dt) = (4.02 ± 1.13) • 10 5 yr. 

It means that the matter streams rather fast from the less massive star to the 
more massive one. In the case of WZ Crv it means that the second component lost 
its mass and the matter streams throw the inner Lagrange point falls into the Roche 
lobe of the second component, which completely consistent with the direct impact 
accretion stream version of a hot spot. 
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Figure 5. O-C diagram. 



VI Conclusions 

• The poorly studied Algol-type eclipsing binary system WZ Crv was investi- 
gated for the first time in two photometric systems V and R. 

• The brightness and color indices were determined at the primary and secondary 
minima and maxima. 

• The brightness at maxima differs, indicating a possible spot in the atmosphere 
of one of the stars. 

• The temperatures, spectral classes and radii of both components were esti- 
mated. 

• The increase of the orbital period was discovered. 
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